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Abstract. We compute the anomalous dimension of the third and fourth moments of the 
flavour non-singlet twist-2 Wilson and transversity operators at three loops in both the MS 
and RI' schemes. To assist with the extraction of estimates of matrix elements computed using 
lattice regularization, the finite parts of the Green's function where the operator is inserted in 
a quark 2-point function are also provided at three loops in both schemes. 
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1 Introduction. 



In calculations relating to deep inelastic scattering the operator product expansion plays an 
important role in allowing one to evaluate the underlying current correlators. In essence there 
are two parts to the formalism. The first is the basis of gauge invariant operators into which the 
current correlators are decomposed and the other is the process dependent Wilson coefficients. 
For high energy experiments the dominant set of operators are those of leading twist defined 
as the difference between dimension and spin. The Wilson coefficients are determined from the 
specific process of interest and are computed in perturbation theory. To understand the physics 
at various energies one requires the solution of the underlying renormalization group equation 
at some particular loop approximation and necessary for this are the anomalous dimensions 
of the operators of the basis. With the increased precision now demanded by experiments 
the goal in recent years has been to ascertain the anomalous dimensions at three loops as an 
analytic function of the moment, n, of the operators. This has been achieved by the magnificent 
computation of |^ |21 El 15 for the twist-2 flavour non-singlet and singlet Wilson operators in the 
MS scheme as well as the Wilson coefficients to the same order which extended the lower loop 
results of 0013 ■ Hence the full two loop renormalization group evolution has been determined. 
However, one feature which cannot be deduced from perturbative techniques is the underlying 
matrix element which is non-perturbative in nature and has to be deduced by use of the lattice 
where various moments of the matrix element have been determined for low moments. See, 
for example, the ongoing work of the QCDSF collaboration, [51 El El El E] > and others, 
|14 | I15 | [THl I17j . However, in making measurements of matrix elements, which ultimately will 
provide accurate predictions for experiments, there are several issues. 

First, to reduce computation time and consequently cost, the matrix elements are deter- 
mined on the lattice in a renormalization scheme geared for this problem which is known as 
the regularization invariant (RI) scheme and its variation called the modified regularization in- 
variant (RI') scheme, |18j . Both are mass dependent renormalization schemes. However, when 
the results are extracted in this scheme they need to be converted to the reference MS scheme 
which is a mass independent scheme. An early example of the application of this approach is 
given in I2()j. Second, to ensure that the results are credible when evolved to large energy 
they must match the perturbative result for the matrix element in the same renormalization 
scheme. There are several ways of achieving this. One is to use a non-perturbative approach 
such as the Schrodinger functional method, 21 . The other is to compute the matrix element 
to as many orders as possible in conventional perturbation theory and then match the lattice 
results to the explicit perturbative results in the same renormalization scheme. Previously that 
has been the approach of Chetyrkin and Retey, 20|, and [22l|2S]. Specifically, various quark 
currents have been considered including the tensor current as well as the second moment of 
the flavour non-singlet Wilson and transversity operators. The latter was originally introduced 
in j24| I25| I2(ij and relates to the probability of finding a quark in a transversely polarized nu- 
cleon polarized parallel to the nucleon versus that of the nucleon in the antiparallel polarization. 
The results of j2()| I22L l2,Sj have been important in the matrix element lattice computations of 

El El El El El Ej • Necessary for the three loop perturbative calculations of the Green's 
functions with the operator inserted has been the full three loop renormalization of QCD in 
the RI' scheme in an arbitrary linear covariant gauge, j201 . Given that there has been an 
advance in computing technology in recent years, it transpires that it is now feasible to measure 
higher moments of the underlying matrix elements on the lattice since essentially there has been 
a significant improvement in numerically isolating a clear signal. Therefore to assist the full 
matching procedure in the ultraviolet region to produce accurate estimates, it is necessary to 
extend the approach of |22| I23j to higher moments of these two classes of operators. This is 
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the purpose of this article where we will consider the third and fourth moments of the flavour 
non-singlet Wilson and transversity operators in an arbitrary linear covariant gauge. There will 
be two parts to this venture. The first is the determination of the anomalous dimensions of the 
operators at three loops in MS and RI'. Whilst the second is to produce to the same loop order 
the value of the Green's function involving the operator itself inserted in a quark 2-point func- 
tion in both schemes. Although it is worth noting that for the lattice, only the Landau gauge 
results are relevant since that is the gauge in which lattice measurements are made. The full 
arbitrary gauge calculation, being more complete, is actually important for internal checks on 
the loop computations. Whilst the even moments are appropriate for deep inelastic scattering 
experiments, the odd moments are accessible on the lattice and can serve the purpose of a test- 
ing ground for technical issues for higher moment lattice matching. Whilst the Wilson operator 
anomalous dimensions are known already at three loops in MS, El) we note that what is 
required are the values of the specific Green's function with the operator inserted which has not 
been determined previously. In also considering the transversity operator, we will deduce new 
anomalous dimensions at three loops in both MS and RI' beyond the earlier two loop calculation 
of [271 1281 QUI I3U1 131j . Moreover, since we will be using symbolic manipulation and computer 
algebra tools and given that the Wilson and transversity operators are both formally similar, 
the actual computations are efficiently performed through the same computer programmes. 

The paper is organised as follows. In section 2 we review the basics of the RI' scheme 
and introduce the properties of the operators we will consider. The full three loop anomalous 
dimensions for these operators are given in both schemes in section 3, whilst the same information 
for the underlying Green's functions are provided in section 4 including the restriction to the 
colour group SU (3) and the Landau gauge. Section 5 records the explicit functions which convert 
all the three loop anomalous dimensions between the MS and RI' schemes with concluding 
remarks given in section 6. Several appendices summarize the projection of the Green's function 
onto the basis of independent Lorentz tensors sharing the same symmetry properties as the 
corresponding operator, as well as the construction of the full operator with these symmetries. 



2 RF scheme. 



In this section we discuss the definition of the RI' scheme and properties of the operators 
we are interested in. First, we recall that the standard renormalization scheme is the MS 
scheme, [HH, where the poles with respect to the regulator are absorbed into the renormalization 
constants. Its widely used extension, MS, which is also a mass independent scheme additionally 
absorbs the finite part ln(47re~''') where 7 is the Euler-Mascheroni constant, [SSI- contrast 
the regularization invariant schemes, JH]) are mass dependent schemes which from the point 
of view of the Lagrangian is the scheme where the quark 2-point function not only is rendered 
finite but also chosen to be unity through the RI definition 
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(2.1) 



where is the quark 2-point function with external momentum p, is the quark wave 

function renormalization constant in the TZ renormalization scheme and ^ is the scale intro- 
duced to ensure that the coupling constant remains dimensionless in d dimensions when using 
dimensional regularization with d = 4 — 2e. However, in practice since taking a derivative is 
(financially) costly on the lattice, one takes a variation of this definition ^2.\\ . to define the RI' 
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scheme which does not involve differentiation, UH]) which is 
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Although this is primarily the key to defining the scheme on the lattice as well as the continuum, 
the full three loop renormalization of QCD has been performed in an arbitrary linear covariant 
gauge and colour group in [22j. Additionally part of the four loop renormalization has been 
performed for SU{Nc) in 20 . However, in [221 the other field 2-point functions were defined 
in an analogous way to ()2.2() . Further, by contrast the 3-point functions were renormalized 
according to the usual MS condition. So that those Green's functions were not constrained to 
be unity. Consequently, the relationship between the variables in both schemes were established 
to three loops and specifically are, [22]) 
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where a = g'^/^lQn'^) in terms of the coupling constant g in the definition of the covariant 
derivative D^, and for the arbitrary linear covariant gauge parameter a, 
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where Tp, Cp and Ca are the usual group theory factors defined by 
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for colour group generators T", C(n) is the Riemann zeta function and the scheme of the variable 
is indicated as a subscript. Clearly only in the Landau gauge, where a = 0, are the variables 
equivalent. Although ultimately we are interested in the Landau gauge for the lattice matching, 
we have chosen to compute in an arbitrary hnear covariant gauge as the extra a. dependence, 
evident in expressions such as (|2.4|) . will provide a central checking tool in the loop computations. 
For instance, in a mass independent renormalization scheme the anomalous dimension of a gauge 
invariant operator is independent of the gauge parameter. Though this is not the case for mass 
dependent schemes such as RI' which will be apparent in our explicit results. 

More specifically the non-singlet operators we will focus on are 



rr\V\...Vn 
QlJ,Ui...Un 



(2.6) 



with n = 3 and 4 where a'^'^ = ^[y^, 7*^] and is antisymmetric in its Lorentz indices. The operator 
iS in both sets denotes the symmetrization of the Lorentz indices {z^i, . . . , and the selection of 
the traceless part according to slightly different criteria for both cases. For the Wilson operators 
this is 

V.^uPw"'^'-"'-"'' = . (2.7) 



Whilst for the transversity operators, [2 
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Therefore, the transversity operator has formally one less traceless condition compared to the 
Wilson operators. 

For the renormalization of an operator in a quark 2-point function, which will be either 
(V'(-p)O{^'"''"(0)Vi(p)) or (V'(-p)O^'''- ''"(0)^(p)), the RI' scheme definition is similar to (g^J, 
jl8l 1191 Oni 1221 123j . However, as the operators we will consider will carry Lorentz indices this 
2-point function will decompose into several invariant amplitudes which may or may not have a 
tree (T) or Born term. For those amplitudes which have a tree term, the RI' scheme definition 
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where (p) is the tree part of {Tp{—p)O{0)tlj{p)) and T is the value of the tree term amplitude 
which may not necessarily be unity. In other words there is no a dependence in T. The 
explicit details of our choice of how to construct the Green's functions in terms of a basis 
of Lorentz tensors satisfying the same symmetry properties as the original operator is given in 
appendix A. This summarizes the procedure we will use to extract the renormalization constants 
of the operators which will then be encoded in the associated anomalous dimensions through 
the respective definitions 
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3 Anomalous dimensions. 



Having described in detail the method of renormahzing in the RI' scheme, we now record the 
exphcit three loop results for the anomalous dimensions in this section. In constructing our 
results we made extensive use of the Mincer package, [311 ES]) written in the symbolic manip- 
ulation language Form, 36 . The Mincer algorithm, determines the divergent and finite 
parts of massless 2-point functions using dimensional regularization in d-dimensions. Therefore, 
it is ideal for the current problem since the Green's functions we are interested in are massless 
quark 2-point functions with the appropriate operator inserted at zero momentum. This is the 
momentum configuration which is measured on the lattice. Moreover, since we are concerned 
with operators which are symmetrized in their Lorentz indices and satisfy various tracelessness 
conditions in addition to being flavour non-singlet operators, there is no possibility of mixing 
into other operators. This is an important observation since ordinarily nullifying the momen- 
tum flow through the operator could lead to the inability to correctly determine the projection 
into the full basis of operators. (For a clear exposition on the deeper perils of operator mixing 
see, for example, [SZIO The fact that each of the operators is multiplicatively renormalizable 
avoids this potential technicality. For our three loop computation we generated the Feynman 
diagrams with the Qgraf package, |3H]- Specifically there are 3 one loop, 37 two loop and 684 
three loop diagrams to be calculated. Though for the operators with no three gluon and two 
quark leg operator insertions the latter total is reduced by 14. Finally, the electronic Qgraf 
output is converted into Form input notation and the Form version of the Mincer algorithm, 
|35j ■ is applied to the 724 Feynman diagrams. The actual Feynman rules for each operator were 
generated automatically in Form. First we constructed the object with the same symmetry and 
traceless conditions as the operators we are interested in. The explicit details for each operator 
are given in appendix B. Then we applied an algorithm which systematically substitutes for the 
covariant derivatives and functionally differentiates with respect to the various fields present to 
produce electronic forms of the 2, 3, 4 and 5-point operator vertex insertions. 

Now we provide our results in MS. For completeness we give those for the two Wilson 
operators and note that we found exact agreement with the results first deduced in [TllSllHl lH^lin] . 
These are 
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where we note that throughout the article when the operator appears exphcitly as a subscript 
on an object, it is regarded as a label and the free indices do not endow the object with ten- 
sor properties. Likewise when we indicate the renormalization scheme on a quantity which is 
evaluated in perturbation theory, that means that the variables in which it is expressed, such 
as a and a, are regarded as the variables in the same scheme. The relationship between the 
variables in either scheme is given in (|2.3() and 1)2. 4|) . For the two transversity operators the MS 
anomalous dimensions have not been given previously and we find that 
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There are several checks on these two expressions. First, as we have computed them in an 
arbitrary covariant gauge their final form must be independent of the gauge parameter in a mass 
independent renormalization scheme, which is apparent in (|3.3j) and (|3.4|) . |32|l41j. Second, part 
of each of the three loop terms has in fact already been determined by the large expansion 
in 122]. There the leading order critical exponent corresponding to the anomalous dimension 
evaluated at the non-trivial d-dimensional fixed point of the QCD /3-function was determined in 
d-dimensions using a method that was originally developed to study the perturbative structure 
of scalar field theories, j421 143j . This critical exponent, |^, encodes all orders information on 
the corresponding renormalization group function at 0(1 /Nf). Therefore, if we formally write 
the leading 0{1/Nf) part of the arbitrary n transversity operator anomalous dimension as, |23j . 
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then the leading order coefficient of the polynomial at three loops is given by 
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where Si{n) = Ya=i 1/^'- Evaluating this for n = 3 and n = 4 reproduces the corresponding 
coefficients in 1)3. 3() and (|3.4() respectively. As a final check we note that we have used the 
method of [ll] to perform the automatic renormalization of Green's functions at high loop 
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order. This entails computing the Green's functions as a function of the bare couphng constant 
and bare gauge parameter. The counterterms are then introduced at the end of the computation 
by rescahng by the known renormahzation constants. Therefore, the remaining divergence is 
absorbed by the renormahzation constant associated with the Green's function. Moreover, 
given the way it has been deduced, the non-simple poles in e are constrained to satisfy a specific 
form depending on the lower order simple poles due to the underlying renormahzation group 
equation. This is a non-trivial checking criterion, especially in the presence of parameters 
such as the gauge parameter and group Casimirs, and it is reassuring to record that all the 
renormahzation constants determined for the above anomalous dimensions precisely satisfied 
this criterion. Implicit in this final check is the fact that the already known two loop anomalous 
dimensions of HSl UHl OOl 1^ are correctly reproduced when the n-dependent results are 
evaluated for n = 3 and n = 4. All these checks therefore give us confidence that not only are 
all our expressions correct but also, for example, that the original Feynman rules were correctly 
generated. 

Having established the MS anomalous dimensions it is then relatively straightforward to 
deduce the anomalous dimensions in the RI' scheme. This is achieved by replacing the renor- 
mahzation constants which scale the bare internal parameters and that of the overall quark 
wave function, by those of the RI' scheme and then impose the RI' scheme definition for the 
operator renormahzation, p. 9)1 . As a check on the resulting renormahzation constants, the non- 
simple poles in e also have to satisfy various constraints emanating from the renormahzation 
group equation, similar to those of the MS scheme. We note, for completeness, that these are 
fulfilled. Therefore, we record the corresponding three loop RI' anomalous dimensions are, in 
four dimensions. 
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for n = 3. Further, for the n = 4 operators we have 
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in four dimensions. Clearly they all satisfy the trivial check that the one loop term is scheme 
independent. Though since the RI' scheme is a mass dependent one, the anomalous dimensions 
will not necessarily be independent of the gauge parameter as is clearly the case above. 

Although we have performed the computation in an arbitrary gauge and colour group, for 
practical purposes it is useful to specify the results for SU{2>). Therefore, the MS transversity 
anomalous dimensions are 
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where Tp = 1/2, Cp = 4/3 and Ca = 3 for SU{3). In addition for the RI' scheme we record 
each of the anomalous dimensions in the Landau gauge since that is the gauge primarily used 
in matching to lattice results. We have 
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for the Wilson operators. Whilst 
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4 Finite parts. 

In this section we record the three loop MS and RI' expressions for the amplitudes of the various 
Green's functions we computed to obtain the previous anomalous dimensions. These are essential 
for lattice matching computations which therefore necessitates their tedious presentation. The 
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specific definitions of the quantities t!q {p) are, as noted before, given in appendix A. It is worth 
pointing out that not all the amplitudes have an a independent leading term. 

First, for the Wilson operator with n = 3, we have 

^(1) MS finite/ x 















p2 = 







+ 



+ 



/367 6^,^, 361 7 , 6^,^,\ ^ ^ 

C(3)a H a + -o? - -C(3) CfCa 

\ 30 5^^ ' 90 9 5^^ )jtA 



+ - 



1087 37 



1 



24 



120 18 



- -« + -a' + -C(3) - -^^TfCf 



and 

y{2) MS finite, n 



+ 



36151 112^,^, 

\ C(3)a 

243 15 ^ ' 

/ 504013 22735 
+ -TTT— a 



9 5 

14671 



25 



810 

224 



a + ^C(3)) NfTpCpCA 



+ 



V 6480 
/ 6480923 

V 19440 

47 



1944 
3727 
' T20" 



4210 



C(3) ) NfTpCl + ^Zl^NfTlCp 



243 



759413 

C(3)a + 2C(5)a + a 



12960 

401 n 

H a 

216 



2563 



C(3) 



67 



+ 



+ 



24 
2811619 



C(5) CpCl 



25 



6480 



+ yC(3)a + 8C(5)a 



51839 
1215 



■a 



C(3)a' 



613 



-a 



162 

/ 28855943 
V 155520 



979 



784 



15 



C(3) + ^C(5) C'pCa 



4C(3)a + 



218971 



539 



a H a" 

15552 162 



11 o 860,, , 
54" 



272C(5) C% 



+ 0(a^ 



(4.1) 



+ 



+ 

/ 86597 



1 . ^ 
1^ + 6"' 



+ ^^(3)« + + - -C(3) ) C^C^ 

V 4860 5^ ' 360 72 5 ''^ M ^ ^ 



+ 



1471891 401 



12 



155520 



—a + —a' + — C(3) C^F 



32363 



+ 



3888 

30365437 



NfTpCp 
68 



216 



36 



209952 

3577 



+ 



+ 



243 
/ 1019471 
V 2099520 " 

4166 



1474 

+ 45^^')"-l0^" 
100 



9 

100 



C(4) NfTpCpCA 



26413 

— C(o)a H a 

27 ^ 2592 

/ 1227463 400 



+ 



52488 243 
/ 208545851 _ 3721 
V 1049760 



C(3) N/T^Cf 



720 C(3)«-^C(4)a + ^C(5)« 
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+ 



+ 
+ 

+ 
+ 



390323 23 2 1 2 

a C(o)a C(4)a 

25920 720^^ ^ 16^^ ' 



17^, , 2 ciuy»^ 



720^^ ^ 16- 
30983 2 475 , 

8640 9^^ ^ 864 



25920 
31163 

972 IH:^ 'i J 

781692217 697 _ 8 _ 57691 

8398080 1a^^ ~ 3^' ~ 1620 

7 , 2 7895 2 ^ s Is 743,, , 

— C(3)a^ - -C(3)a^ H C(3) 

18^^ ' 2592 3^^ ' 24 540^^ ^ 



■C(3) + ^C(4) + ^C(5)) CfCI 
2217 697^,^, 8^,^, 57691 



a 



f C(4) - ^C(5)) C'^Ca 



1161367 733,, , 458737 25 , 2 

C(3)a a C(3)q;^ 

43740 54 ^ 20736 9 ^ 
«T7 o Q1 noriQ 



54 

617 2 2 , , o 
|C(4) + 24C(5)' 



31 o 9203 



216 



972 



o3 + 0{a^ 



(4.2) 



For the n = 4 case we have 

^(1) MS finite / X 



^ /1871 4 ,^ 



+ 



+ 



/26869109 
V 324000 



+ 



345682991 43553 



6480000 
6041063 



NfTpCp 



162000 
32796795659 
43740000 



3600 

2 



a + |o^ + fc(3))c^ 



a 



46,,„, 48917 

H C(3)a a 

3 1296 

628 



+ 



+ 



+ 



2025 ' 15 
/ 63233459093 628 



C(4) NfTpCpCA 



V 437400000 
77018 



2010352 

— C(3)a-\ a 

45^^ ^ 30375 



628 

+ -67^^^') + T^^(') ) "^^^^^^ 

/ 1335574847 2512 

+ 



V 10935000 ' 405 
/ 1578785326351 43861 



C(3) ) iv/r|CF 



1399680000 
56000717 



720 C(3)a-^C(4)a + ^C(5)« 



+ 



414720 



-a 



1753 
360 



C(3)q2 



16 



C(4)a' 



3^^ ^ 34560 



986237 2 1.. X s .00^ o 



7859 



-a 



+ 



2350679 ,, , 
1552257600373 



1728 
((5)) CfCX 



1749600000 

181148459 
777600 



1200 ' 3 
+ 1555£«3)« + 4a5). 



2250 

32 

a + — C(3)a' 
15 



1336787 
51840 



-a 



12 



+ 



1739 



205 
192 



a" + 



395693 
2700 



C(3) 



50 ^^^^ + ^^^^^ ^ ^^^^ 



/ 203923883969 
V 2916000000 



29329,, , 685201111 
C(3)q H a 



450 
911003 



157, 2 
T^^^^)" + M800 



4860000 

565 

+ -C(3)a^ - —a" 



2 
3' 



864 



40500 



75 



(4.3) 



and 

^(2) MS finite , , 











{—- 


h —a I 




V160 


32 y 



+ 



+ 



/ 70559 3 , 4991 

C(3)a H a 

V 115200 40^^ ^ 23040 



+ 



40907 7453 



1 



96000 



57600 

2 



a 



384° +5«^"'^^ 



+ 
+ 



119 

^N,T,C, a- 

864000 120^^ ' 

/ 28776619 4054403 77,, 



104179 53 , ^ ^rrr. ^ 

503680° +2i6<P))WFC. 



V 8640000 
/818118521 



4003 5 , 22103567 
2i6o^(3)« + ^C(5)a + -gg3^a 



+ - 



V 55296000 

6641 , 2 1 2 69635 o 913 3 

- C(3)a H C(5)a H H a'^ 

34560^^^ 96^^^ 110592 9216 

198233 263 \ 2 

" ^8800"^^^^ + -96"^^^^ ^^^ 

318424073 7129 , 1 , , 20896549 

\ C(3)a + -C(5)a a 

25920000 12000^ ' 2^^ ' 6912000 



80^^ ' 



502243 



29 



-a 



1024 

^2 



+ 



1382400 

- mh))ClCA + ^^NfnCF 
^ ' " ^ 54000 ^ 

/ 37693459 13 .,0^, 239928131 
V 345600000 ~ 60^' ^ 207360000 



1396697 
108000 



C(3) 



a + 



27433 
115200" 



-a 



41 o 272419 , , 38 , , , 

a C 3) + — C(5 Ci 

4608 27000 ''^ ^ ^ 3 ''V M F 



Turning to the case of the transversity operator, for n 

1 

= ;/.2 18 



(4.4) 



3 we have 



y{l) MS finite/ x 



1 + 



109 5 



+ 



48941 
810 
/ 26467 
V 810 



rC(3)a 



1223 



360 



a 



+ f C(3) ) CfCa 



119 



17 



+ -T^" - - ^C(3) Cp + 



18 



36 



2197 



NfTpCp 



13 



+ 



/2254181 124^, , 32359 

C(3)a-\ a 

V 4374 15 ^ 1620 

+ ^C(3) + ^C(4))W.C. 



+ - 



+ 



+ 



81 

699751 



104^,^, 18349 



+ - 



C(3)- 
416 



C(4) NfTpCi^ 



C(3) ) AT/rlCp 



11176. 104 
135 
177970 
2187 

260680459 5537^,^, 3^,,, 3^, , 

\ C(3)a + -C(4)a - -C(5)a 

349920 180 ^ 8^^ ^ 2^^ ^ 

95953 133^,^, o 3 o 4^, , , 

1296 80 ^ 16^^ ^ 3^^ ^ 
35543 9 1,, , o 2227 o 

2160 3^ ^ 864 
194707^,^, 463^,,, 139^,^,\^ ^ 



1296 
/ 77900401 
V 174960 
5 
6 



48 



6 



3983^,^, , 635981 
-C(3)a + 4C(5)q + a 



90 

K(3)«^ + ^"^ + C(3)a^-:^«' 



4860 
1 

2"4' 



9163 
T35" 



C(3) 



+ 22C(4) - 20C(5)) C'fCa 



/265849 410 74771 ^26 2 

C(3)a a-\ C(3)a 

V 7290 9 ^ ' 972 3 ^ ' 



5063 2 2^,^, 3 115 3 
- 3^^^^" + 216" 



15518 



648 

f C(4) - 32C(5)') C| 



405 



C(3) 



(4.5) 



with further calculation giving 
y{2) MS finite, N 



^(3) MS finite. ^ 



y{l) MS finite/ x 



3 (1) MS finite, . 
2 tpaP-" D" DPtj) 



(4.6) 



For the fourth moment of the transversity operator we have 

293 



y(X) MS finite / s 



1 

32 



1152 384 



— -a Cpa 



+ 



+ 



/ 2037217 3127 397 2 

1 a H a 

V 829440 18432 9216 

_^ _ 8099 1595 29 



13 
32' 



C(3) CfCa 



5184 
118621 



4608 



a + 



4608 



rC(3) 



NfTpCp 



103680 
2440299949 



a 
59 



H C(3)a 

111974400 160^^ ^ 



384991 
414720 



a 



14681 
6480 



C(3) - -C(4) j NfTpCpCA 



14 



+ 

+ 
+ 



/ 421 19449 4^,„, 2369273 

C(S)a-\ a 

V 11197440 9^^ ^ 1244160 

+ ^C(3) + tc(4)) NfT^Cl 



1080 
/ 25433893 



+ 



16 



3583180800 
89697707 

26542080 " 
17 



+ ^C(3) ] Nf^nCp 
C(3)q - 



V 6998400 ' 81 

112663077941 99239 



69120 
2879 



— C(4)q + — C(5)a 
256^^ ' 96^^ ^ 



^, , 2 541433 

C(5)a H 

384^^ ^ 737280 

503123 



2 1 . ^ 12979 

H C(3)a H 

96^^ ^ 110592 



a 



+ - 



82944 
2225988241 
89579520 



10093^, , 



16376411 
2488320 



-a 



360 
/ 7904501 



477515 2 1 X s 
663552" -32^^')" " 

|C(4) + |C(5)) C^C^ 



323 
12288 



a 



295^, , 1045465 
144^^^)" + ^48832"" 



1 



V 3732480 

64399 2 , 1 y-ro\ 3 



rC(3)a2 



+ 



165888 48^^ ^ 



1007 



55296 



13739 
12960 



C(3) 



(4.7) 



The remaining amplitudes are related to the first through 

4 



y{2) MS finite / s 
y{Z) MS finite / ^ 



_ _ ^ (1) MS finite / X 
_ 1 ^(1) MS finite / ^ 



(4.8) 



Finally, for practical purposes we provide the general results for the specific case of the 
Landau gauge when the colour group is SU{3). Therefore, we have 



^(1) MS finite. ^ 



5;7(3),a=0 



4 

9" + 



50 4432 56 ,, / 



+ 



4210 



729 



iV/+(- 



1665047 



+ 



279011797 
131220 



7290 
1717789 
2430 



- f C(3)) N, 



<(3) + ^«5) 



+ 0(a4 



(4.9) 



and 



y{2) MS finite/ x 



5C/(3),a=0 



1 

-3 + 



214 



+ 



81 

/ 1227463 



a + 



32363 4763093 



5832 



87480 



152 



C(3) 



157464 729 



+ ^(3))^, 



V 



15 



+ 
+ 



1364405723 
4723920 

/ 8619089351 



1000 



C(4) Nf 



V 4723920 
8599 



12125507 
32805 



C(3) 



+ 



972 



C(4) + ^C(5) 



+ O(a^) . (4.10) 



For the RI' scheme we note that 



y{l) RI' finite/ n 



5(7(3), a=0 



4 

9° + 



50,, 44168 56,, , 

—Nf + + — C 3) 

27 ^ 1215 15^ ' 



+ 



4210 



729 



2738978 28 



10935 



C(3) Nf 



+ 



326345791 
131220 



2430 



27 



y{2) RI' finite/ s 



5i7(3),a=0 



(4.11) 



Further, for the next Wilson operator 

SU(3),a=Q 



y(l) MS finite / N 



and 



y{2) MS finite / n 



SU{i),a=G 



As yf^^ ^""^ (n\ 

■4)'ylJ'D''D'^DP'4)^^' 



SC/(3),a=0 



7484 

- 1 H a + 

675 

/ 1335574847 

+ 



6041063 ,, 431713457 524, 

Nf + C 3) 

243000 ^ 1822500 15 ^ ^ ^ 



2512 



V 32805000 1215 
1349;]8888G469 

+ ■ 
+ 



C(3) ) Nf 



984150000 
706189399771421 



1215 
112503104 



27 



78732000000 

43507^,^, 7180^,^, 



C(3) 

54675 ' 
+ 0{a*) 



(4.12) 



119 731129 23 

1200 ^ ^ 432000 ^ 90^^ ^ 



51959 
162000' 



N, 



f 



232632277 



19440000 



ic(3))A'; 



+ 



1047728166241 
9331200000 



109467991^^^^ _^ 13111 



2916000 



648 



C(5) 



+ O(a^) . 
(—1) by construction, we note 



(4.13) 



y{2) RI' finite / x 



5(7(3), 0=0 



40 



+ 



a + 



51959 
162000 



119 850873 23,, , 

Nf + + — C(3 

1200 ^ 432000 90^ ^ 



Nf + 



266088707 
19440000 



ic(3))iV, 



+ 



+ 0(a^ 



435587120587 
3110400000 



— C(3, + Hllic,5) 



583200 



648 



a 

(4.14) 
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For the transversity cases, when n = 3 we have 



^(1) MS finite, ^ 



5C/(3),a=0 



1 

18 



218 



■a + 



4394 



'-N, 



+ 



243 

177970 416 
6561 ~ 243 



669202 
3645 

2 



15 ^ ' 



C(3) iV/ 



/98639141 12712^., 1040^,,,,,, 



V 98415 
1020141085 



1215 
59050063 



27 



and 



^(1) RI' finite, X 



324 

SC/(3),a=0 



157464 ^ 43740 



C(3) 



ijjaP'^ D" DPi) 

Finally, for the transversity moment n = 4 we have 



27 



18 ^ ' 



+ 0{a^) (4.15) 



(4.16) 



^(1) MS finite 
il)cri^'^D'^DPD^^^^^ 



SU(3),a=0 



1 293 

— + a + 

32 864 



118621 13151593 85^, / 



+ 



155520 ' 1866240 72 

/25433893 16 A , . 

+ C 3) N? 

V 20995200 243^ V ^ 

, ,____l?l!c(3)-|c(4))iV, 



503884800 



1944 



+ 



2013899793847 146176079 



8062156800 



2799360 



C(3) 



2693^,,, 52991 



with clearly 



^(1) RI' finite . . 

ipaf^-'D'^DPD^tp^^' 



3456 

SU{3),a=0 



2592 



a3 + 0{a' 



(4.17) 



(4.18) 



5 Conversion functions. 



An additional check on our computations is provided by the conversion functions for each of the 
operators we have considered. These functions allow one to convert the anomalous dimension 
of the operator in one renormalization scheme to that in another scheme and are defined by the 
ratio of the renormalization constants in both schemes 



Then, IS?], 



To («RI') 



7RI' 

Co(a.a) = ^. 



To ^ («Ms) - ^ ("Ms) Q— 1" («MS' «MS 



MS 



- «MsT^^ («Ms) ( 



""MS'^MS 



(5.1) 



(5.2) 
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where one needs to express the MS variables in terms of the RI' scheme using (|'2.;-i|) and (|2.4j) 
in order to compare with the anomalous dimensions from the explicit computation. We record 
that the conversion functions for the various operators we are interested in here are 

CpCL 



a 



(a, a) 



1 + (27a + 107) 



18 



+ 



(864000^ - 93312C(3)a + 409104a - 715392C(3) + 3302464j Ca 
+ (60480a2 + 327600a + 373248C(3) + 327549) Cp 



1475960rpiV/] 



+ 



G 



(a, a) 



51840 

(ll032200a3 - 13915152C(3)a2 - 466560C(5)a2 + 64538856a2 

- 141379344C(3)a + 8398080C(5)a + 319887792a 

- 635381280C(3) + 25660800C(4) + 142767360C(5) 
+ 2356357048) Cl 



+ 



+ 



(4607280a^ + 5785344C(3)a2 + 32256792a^ - 13841280C(3)a 
+ 33592320C(5)a + 180233856a - 297743040C(3) 

- 76982400C(4) - 59719680C(5) + 1051633031) CaCf 
(35085312C(3)a - 97555104a - 75098880C(3) 

- 93312000C(4) - 1625432200) CATpNf 

(2177280a^ - 23328000C(3)a2 + 27799200a2 - 73685376C(3)a 
+ 90339057a + 319139136C(3) + 51321600C(4) 
+ 201553920C(5) - 607345686) C| 

(31104000C(3)a + 63327960a - 303948288C(3) 

- 93312000C(4) + 922104436) CpTrNf 



(13824000C(3) + 250653280) T^Nf 



2799360 



+ O(a^) (5.3) 



1 + (525a + 1871) 



CpCL 



+ 



(l8315000a^ - 23328000C(3)a + 88528500a - 103680000C(3) 

+ 603802180) Ca 
+ (21330000a2 + 119182200a + 62208000C(3) + 98349057) Cp 

-264322520T.A',l^ 

(239334750000a^ - 340977600000C(3)a2 - 2916000000C(5)a2 

+ 1428857212500a2 - 3356364600000C(3)a 

+ 174960000000C(5)a + 7142849746875a 

- 12700608624000C(3) + 335689920000C(4) 

+ 2504844000000C(5) + 48069511158775) Ca 
+ (229047750000a3 



142300800000C(3)a2 + 1689792435000a^ 



- 1524733632000C(3)a + 839808000000C(5) 
+ 9165889557000a - 1770530400000C(3) 



a 
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- 1007069760000C(4) + 653184000000C(5) 
+ 18744964493980) C^Cf 

+ (816480000000C(3)a - 2158137000000a - 1801163520000C(3) 

- 1464998400000C(4) - 31372620527200) CATpNf 

+ (l45435500000a3 - 366249600000C(3)a2 + 137261 14200000^ 

- 1048904640000C(3)a + 3148063990200a 
+ 4800009888000C(3) + 671379840000C(4) 
+ 3359232000000C(5) - 8872064364708) C|, 

- (488332800000C(3)a + 2684380392000a - 4552485120000C(3) 

- 1464998400000C(4) + 18121905428720) CFTpNf 

+ (217036800000C(3) + 4952079510400) T^N?'^ 



34992000000 



+ 



+ O(a^) (5.4) 
1 + (33a + 109)%^ 

io 

+ [(6600a2 - 7776C(3)a + 32067a - 47952C(3) + 228974) Ca 
+ (6480a^ + 30900a + 31104C(3) + 10917) Cp 

(3407670a^ - 4575204C(3)a2 - 58320C(5)a^ + 20121777a2 

- 46022256C(3)a + 2799360C(5)a + 100170405a 

- 196675020C(3) + 3732480C(4) + 45081360C(5) 
+ 693358478) Cl 

+ (2334420a^ - 46656C(3)a2 + 159563520^ - 12324960C(3)a 
+ 11197440C(5)a + 86296752a - 34434720C(3) 

- 11197440C(4) + 214883180) CaCf 
+ (11384064C(3)a - 30726648a - 17280000C(3) 

- 24261120C(4) - 463372640) CaTfNj 

+ (l399680a^ - 6065280C(3)a2 + 13024800a2 - 13965696C(3)a 
+ 26888652a + 108183168C(3) + 7464960C(4) 
+ 22394880C(5) - 153974772) 
- (8087040C(3)a + 27287280a - 69133824C(3) 

- 24261120C(4) + 231549328) CpTrNf 

+ (3594240C(3) + 70515200) T^N? ^^^^ + 0{a'^) (5.5) 



699840 



1 + (75« + 293)^ 
36 

^ r / _ . _ _ _ o 



(64890a2 - 77760C(3)a + 309195a - 414720C(3) + 2302897j Ca 
+ (63720a^ + 380940a + 207360C(3) + 40494o) Cf 
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+ 



- losgessTpiV/] 

( 6771276000^ 



25920 

918993600C(3)a2 - 18662400C(5)a2 

+ 4008299580a2 - 9063653760C(3)a + 466560000C(5)a 
+ 19846116045a - 36380232000C(3) + 783820800C(4) 
+ 8939289600C(5) + 140182687541) 
(517071600a^ - 102643200C(3)q2 + 38406474000^ 

- 3332482560C(3)a + 2239488000C(5)a + 21797212320q 

- 9369146880C(3) - 2351462400C(4) 

+ 447897600C(5) + 58907275400) CaCf 
(2217093120C(3)a - 6005931840a - 6177116160C(3) 

- 4777574400C(4) - 94522187168) CaTeNj 

(279936000a^ - 1194393600C(3)a^ + 30138480000^ 

- 4503859200C(3)a + 8684647200a + 18380113920C(3) 
+ 1567641600C(4) + 2985984000C(5) - 24416900640) C| 

(1592524800C(3)a + 6750907200a - 16028098560C(3) 

- 4777574400C(4) + 57917250880) CpTpNf 



+ (707788800C(3) 
+ O(a^) 



15192521216) T|iV/ 



111974400 



(5.6) 



where the couphng constant and gauge parameter are in the MS scheme. We note that the same 
RI' anomalous dimensions are determined as previously. 



6 Discussion. 



We have provided the finite parts of various Green's functions required for the renormahzation 
of the n = 3 and 4 moments of the non-singlet twist-2 Wilson and transversity operators at three 
loops in both the MS and RI' schemes. Since these are available at several loop orders, the hope 
is that they will be central to the extraction of accurate values for the matrix elements which 
will be measured on the lattice. From another point of view the new MS anomalous dimensions 
which are now available for the moments up to and including 4 for the transversity operator will 
provide a useful check on the full n-dependent three loop transversity anomalous dimensions 
when they are eventually determined. The impressive symbolic manipulation machinery which 
achieved the arbitrary n anomalous dimensions for the twist-2 flavour non-singlet and singlet 
Wilson operators, m^lElE]) can be applied to the transversity case. Whilst this can be achieved 
in principle, in the interim one could follow a similar direction to the earlier approach of |39| I40j 
where the anomalous dimensions of the Wilson operators were determined to moment n = 10 
and later to higher moments, n < 16 (except n = 14), |45| I46j . These explicit moments were then 
used to construct solid approximations to the full anomalous dimensions which were then shown 
to be credible in a substantial range of the x variable in, for example, |32I- Given the advance 
in computer capabilities since ^2 , 40 , it would seem to us that a fixed moment computation 
for the anomalous dimensions of the higher moments of the transversity operator is certainly 
viable. Moreover, one would not be constrained by the choice of an arbitrary covariant gauge 
and the use of the Feynman gauge would therefore reduce computer time. 
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A Projection tensors. 



In this appendix we record the decomposition of the various Green's functions we have computed 
into the various tensor bases. For the Wilson operators there are two independent tensors 
consistent with the symmetries and tracelessness conditions of the original operator. By contrast 
the Green's functions involving the transversity operator require three independent tensors. 
First, for the Wilson operators, for n = 3, we define 



(V'(-p)or(o)V^b)) 



{d + 2) 

2 



(p2)3 



{d + 2) 



{d + 2) 



1 



{p- 



,2^3 



(A.1) 



Then 



{d + 2){d + A) 
4(d2 - 1) 

_ (d + 2y 
4(d2 - 1) 
1 



p^'p^'p'^i) - 



{d + 2) 
2 



?7^- p 



{d + 2) 



p- 



{d + 2) 



and 



X j^{i^{-p)o^''{omp)) 

{d + 2 



(A.2) 



+ 



4(d2 - 1) 

{d + 2)p^ 
12{d? - 1) 



p^p^'p'^i - 



p^ 



{d + 2) 
2 



{d + 2) 



{d + 2) 



(p2) 



,2U 



{^P{-p)O^''{0)i}{p)) . 



(A.3) 



In this appendix and the next in order to compactify our expressions, we note that the sym- 
metrization on Lorentz indices using parentheses excludes the standard normalization factor. 
So, for instance, 



(A.4) 



Further, in our convention when this operation involves a tensor which is itself symmetric, then 
one only counts the independent object once. So, for example. 



(A.5) 



For the n = 4 moment of the Wilson operator we decompose the Green's function into two 
independent tensors with 

(^(-P)or^(o)v;(p)) = E^i 



{p)wr'{p) 



(A.6) 
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Then each individual amphtude is given by the projections 

1 



4(ci + 4)(d + 2)(d2 - 1 

(o)V(p)) 



4(d + 4)(d2 -1) 

fj,uap mip)) 



4(d + 4)(d2 - l)p2 

(d + 3) 



4(d + 4)(d2 - l)p2 



(A.7) 



where the two basis tensors are 



2p- 



2 p'^'-p p T 



(d + 2)(d + 4) 



wr'ip) = 



(A. 



Finally, for this particular case to ensure that the tree term of one of the amplitudes involves 
unity, we have chosen to redefine the amplitudes as 



tjj'yt^D'^D''DP4> 
y>(2) 



ip) 



(d + 2)(d + 4) 



(1) 



ip) 



1 



-s 



(2) 



2(d + 4) i^i^^ D'^ DPijj 



{P) 



.(2) 

''tp-yP'D'^D'^DPrl} 



ip)- 



(A.9) 



Turning to the cases for the transversity operators then for n = 3 the decomposition is 

3 



m-p)or'mm) = j:^.^D<^DP^m 



=1 



ipaP-'^D'^DPip^ 



ip) ■ 



(A.10) 



Hence the coefficients are given by the three projections of the Green's function 



.(1) 

''il)aP-''D''DPil) 



^(2) 



ip) 



12(^2 -4) (d2 - 1)^1 

1 



■^{S)puap 



iP) - 



-,(3)puap 



12{cP - 4){cP - ly 



ip) 



12(d2-4)(d2-l) 1 

jd^ + £-2d + 4) iZ)^uap 
12d((i + 4)(d2-4)((i2_l) 2 



ip) 



-,{3)pi/ap 



.(3) 



ip) = 



12d(d + 4)(d- 1) ^ 
1 



ip) 



puap mip)) 

jd^ + 2) 
12d(d + 4)((i-2)(d2 -1 



12d(d + 4)(d 

X (V(-;>)CTM.ap(o)v;(p)) 

where the choice of basis tensors is 

.(3W,^, _ _ id + 2)id + 4) ^xppupa^pp^ ^ ^ 2)a^^r?('^V^PA 



(A.ll) 



(p) = - 
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/^(f crp) 2 



{3)fj.uap 



(P) 



Finally, for the n = 4 moment of the transversity operator we set 

3 



{ij{-p)Ol^^''''\Q)m) 
where the three amplitudes are given by 



i=l 



(A.12) 



(A.13) 
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16{d + 6){d + 4){d + 3){d + l)(d2 - - 2) ^ 
in terms of the three basis tensors 
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(A.15) 



Finally, it is worth observing that in the two loop computation of the anomalous dimensions of 
ini El HH UHl HH Oni the Feynman rules were constructed by introducing a null vector 
which projected out a part of the Green's function which had a non-zero tree term and therefore 
allowed for the extraction of the anomalous dimension. Since the lattice matching specifically 
requires information on the full Lorentz structure of the Green's function we cannot follow that 
approach as the null vector would exclude access to several of the amplitudes which are extracted 
from the data by choosing different momentum configurations. 
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B Operators. 



In this section we list the full form of the four operators we considered which satisfy the sym- 
metrization and tracelessness properties. For the Wilson operators we have 

= lyM- - —L_r^(t^'^w^)\ (B.l) 

for n = 3 where 

Wi^"" = ^il;j^i'D''D''^ (B.2) 

and for n = 4 
where 

W'""'P = ^i[jj^''D''DPD''^ . (B.4) 

For the transversity operators, which have one fewer traceless conditions compared to the Wilson 
operators, the full operators are more involved. For n = 3 we have 

^^IJ-ucrp _ rpixvap _ (c?^ + 4(j + 2) ij,{y\B _ 2 rp{v\lJ-9 crp) 

" d{d + 2){d + Ay <^\^ d{d + 2){d + A) ^l'' 

^ d(d + 4) + d(d + 4) l^l'^ 

d(d + 4)^ l^'^ + (d + 2)(d + 4)^ ''^'^ 

where 

rppuap ^ ^^a^^^D^Df^ . (B.6) 
Finally, for n = 4 the fully symmetrized and traceless operator is 

f^HVapX _ rpp,V(Jp\ _ + 7d + 8) rpp.{y(y\d pX) 

^ " {d+i){d+m+^) ''^ 

_| + ^) rppB Ay<T pX) _ ^ rj.{y\p\a\e pX) 
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where 



yM^apA ^ ^ipaf'^^DPO'^D^'^i; . (B.8) 
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